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Abstract 

In this paper, we obtain a general expression for the entries of the 
/th {I G Z) powers of even order {2k + 1)—diagonal Toeplitz matrices. 
Additionally, we have the complex factorizations of Fibonacci 
polynomials. 


1 Introduction 

Tridiagonal matrices are used in solving problems in a large variety of disci¬ 
plines including mathematics, physics and others. In [T] it indicates that using 
tridiagonal matrices to finding the magnetization dynamics of outline an ap¬ 
proach. Ahmed et al. [2] enquired some explicit formulas for the powers of 
a square matrix. Otele§ and Akbulak [3] gained a general expression for the 
entries of the power of complex tridiagonal matrix, in terms of the Chebyshev 
polynomials of the first kind. A general expression for the power of complex 
persymmetric antitridiagonal Hankel matrices are studied in [4] . Silva [5] stated 
integer powers of real anti-tridiagonal matrices using Chebyshev polynomials. 
The expressions for powers of complex skew circulant matrices are found in [B] . 
Burn and Bozkurt [7] are examined powers of tridiagonal matrices. The powers 
of Toeplitz matrices are calculated in [8] . 

Note that our results in this paper are the more general form of the results 
obtained in [Hilo]. Salkuyeh [9] are calculated powers of the tridiagonal toeplitz 
matrices as /c = 1. One can easily see this in cni, taking fc = 2 in Theorem 1. 

This paper is organized as follows: in the next section, we give the eigenvalues 
and eigenvectors of {2k + 1)— diagonal n— square Toeplitz matrix. In Section 
3, the 1—th power of {2k -f 1)— diagonal Toeplitz matrix we will get by using 
the expression {AkY = Pk{JkY{Pk)~^ [12], where Jk is the Jordan’s form of Ak, 
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Pk is the transforming matrix. In Section 4, numerical examples are given. In 
Section 5, determinant of (2k + 1)— diagonal Toeplitz matrix obtain complex 
factorizations for Fibonacci polynomials. 

Let Ak be the (2k + 1)—diagonal n—square Toeplitz matrix as following 
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0 ■ 
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where a S R and 6, c S R. \ {0}. 

For example, is given as follows: 


^5 = 


a 0 
0 a 
0 0 
0 0 
0 0 
c 0 
0 c 
0 0 
0 0 
0 0 


0 0 
0 0 
a 0 
0 a 
0 0 
0 0 
0 0 
c 0 
0 c 
0 0 


0 b 
0 0 
0 0 
0 0 
a 0 
0 a 
0 0 
0 0 
0 0 
c 0 


0 0 
b 0 
0 b 
0 0 
0 0 
0 0 
a 0 
0 a 
0 0 
0 0 


0 0 
0 0 
0 0 
b 0 
0 b 
0 0 
0 0 
0 0 
a 0 
0 a 


( 1 ) 


2 Main Results 


Theorem 1 Let A^ he n—square (2k + 1) —diagonal Toeplitz matrix as in ([ip. 
Then the eigenvalues and eigenvectors of the Ak are 

X^ = a-2V^ccos(^),r=T^ ( 2 ) 

\n + k J k 

for n = 2ks (s S N) and 


{Pk)j = 


(fc-i) 

Uo(ar) 0...0 


(fe-l) (/c-l) 

0...0 fj,^^ U n-k (ar) 0...0 


(3) 
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j = 1,1 + fc, 1 + 2fc,..., n - fc + 1; r = ^ 


(Pkh = 


(fc-l) (fc-l) (fe-2) 

0 Uoiar) o”^ fi^Uiiar) ... o”^ fi^Un-k{ar) o”^ 

^ (4) 

-2 
k 


j = 2,2 + k,2 + 2k,...,n-k + 2-r = 


(/c-l) (fc-1) (fc-3) 

iPk)j= Q Q Uo{ar) o”^ lJ.^Ui{ar) ... 0^?^ /i^C/n-fc(ar) o”^ 

^ (5) 

j = 3,3 + fc, 3 + 2k,... ,n — k + 3;r = ^ 


(Pkh = 


(k-2) (fe-1) (fe-1) 

0...0 ?7o(ar) 0...0 ^^Ui{ar) ... 0...0 Un-k (ar) 0 

k 

( 6 ) 


n T 


J+1 


j = fc — 1, 2A: — 1, 3fc — 1,..., n — 1; r = ^ 


i.Pk)j = 


(k-l) (fc-1) (fc-1) 

0...0 Uo{ar) 0...0 fj,iUi{ar) ... 0...0 U n-k (ar) 


j = k,2k,3k,... ,n-,r = | 

Ar —a 


(7) 


where fi = ^,ar = and Un{.) is the nth degree Chebyshev polynomial of 

the second kind. 


Proof. Let 


F.. := 


X — a b 
c X — a b 

c X — a b 


c X — a b 
c X — a 


For initial conditions det (Fg) = 1 and det (Fi) = X — a, we have 
det (F„) = (A - a) det (F„_i) - 6cdet (F„_ 2 ). 
The solution of difference equation in ([5]) is 

det(F„) = (6c)^C7„(0) 


( 8 ) 


(9) 


( 10 ) 
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where 0 = 
kind [ 11 ]: 


and Un (■) is the nth degree Chebyshev polynomial of the second 


Un (x) 


sin{{n + l)arccosx) 
sin{arccosx) 


All the roots of {7„ {x) are included in the interval [—1,1]. Let 


|A/ra — Afcj — (A) 


( 11 ) 


and 



r det(i^.)2-det(An_i)4 

if k = l 


II 

<1 

det(Ft)"-(*-i)'= det(Ft_i)‘'=-”, 

1 if n < tk, 1 < k < 

^ and t gZ 


[ det(T2)"-'=det(Fi)2'=-", 

if ^ < k < n — 1. 

( 12 ) 

So, the equality (IT^ is written as 




r ( 6 c)t[C/|( 0 )-[/|_^( 0 )], 

if k = l 


II 

<1 

(6c)t[/r^‘-'^"(0)4‘^r(0), 

if n < tk, 1 < k < ^ 

and t G Z 


[ {bc)^U^-'^{9)U^^-^{9), 

if ^ < k < n — 1 



(13) 

where 9 = ^^ 7 =. There are the relations between 14 (.), 114 (.) and 17„(.) 
polynomials as following 


v^{e) =u^{9)-u^.,{e) 

W^i9) =t7n(0)+[/n_i(0) 


(14) 


here 14 (.) and 114 (•) are the nth degree Chebyshev polynomial of the third 
and fourth kind, respectively m- Substituting m into m, we possess 


r {bc)^v^ie)w^{e), 

Aa. (A) = \ {bc)iur^*-^^''i0)u*-~r{9), 

[ (6c)t[/r"(^?)c/i^"-"(0), 


The eigenvalues of Ak obtained as 


if k = l 

if n < tk, 1 < k < ^ and t G Z 
if f < fc < n — 1 . 

(15) 


Xr = a — 2'/bc I 


krn 


,forr = 1 , 


The multiplicity of all the eigenvalues A^ (r = 1,2,..., -1) of the matrix Ak are 
k. Since rank{XrIn — Ak) = n — k, for each eigenvalue Xr correspond ^ Jordan 
cells Jr{Xr) in the matrix Jk- That is 


n/k 


njk 


Jk — diag{Xi ,..., Ai, A 2 ,..., A 2 , 


njk 



(16) 
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Consider the relations {Pk)~^AkPk = Jk [12] we have to obtain the matrices 
Pk and {Pk)~^ and derive the expression of the matrix {AkY for I G N. Let us 
denote j-th column of Pk by {Pk)j {j = 1,..., n). Then 

AkPk = ((Pfc)iAi ... (Pfc)f Ai ... (Pfe)(^A^ ... (Pfc)„A^). (17) 

From Eq. (EZD, we have the system of linear equations as follows: 


Ak{Pk)i 

(Pfe)iAi 

Ak{Pk)^ 

= (Pfc)fAi 

^fe(Pfe)^+i 

= (Pfe)|. + lA2 

AkiPk)s^ 

= {Pk ) 2" A 2 

k 

Ak(Pk) (r,-l)n 

- T: - 

= (Pfc)(n-l)n As 

Ak{Pk)n 

= {Pk)n^^ 


Solving the set of the system of linear equations in (fTSl) as regards 


iPk)j = 


(fc-i) (fc-i) (fc-i) 

Uoiar) O”^ f^^Uiiar) ... o”^ o”^ 


1 T 


j = 1,1 + fc, 1 + 2fc,..., n - fc + 1; r = ^ 


(Pnh = 


(fc-1) 


(k-1) 


(k-2) 


0 UQ{ar) 0...0 fi^Ui{ar) ... 0...0 {ar) 0...0 


j = 2, 2 + fc, 2 + 2fc,..., n — fc + 2; r = ^ 


(Pnh = 


(fc-1) 


(fc-1) 


(fe-3) 


0 0 Uo{ar) 0...0 fi^Uiiar) ... 0...0 fi'^Un^^ar) 0...0 


j = 3, 3 + fc, 3 + 2k,... ,n — k + 3;r = ^ 


(Pnh = 


{k-2) (fe-1) 

o”^ Uoiar) 0^?^ ^i^Ul{ar) 


(fc-1) 

0...0 fj,^^ U n-k (ar) 0 
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j = fc — 1,2A; — 1, 3fc — 1,..., n — 1; r = 


{P^)3 = 

j = k,2k,3k,...,n;r = j: 

where = ^,Oir = and Uni-) is the nth degree Chebyshev polynomial of 

the second kind. ■ 


(fc-i) (fc-i) (fc-i) 

0...0 Uoiar) 0...0 fi^Uiiar) ... 0...0 fi^^Un-k (ar) 


3 The integer powers of the matrix Af. 

Considering ®, ®, ® and ®, we write down the matrix 


Pk = 


Uoiai) 

0 


0 

AUiiai) 

0 


0 0 

Uoiai) 0 


0 

0 

AUiiai) 


U n-fc fai ) 0 0 

0 fi^^Un-kiai) 0 


Uoia2) 

0 


Uoiai) 

H^Uiia2) 

0 

AUiiai) 

/r^C/n^(a 2 ) 

0 

U „ fc (a-f ) 
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C/o(a^) 

0 


0 

C/o(a») 


0 

0 


0 

0 


^i2Ui{a^) 


0 


^ 17 ^ U „ k (an) 


Uo{a^) 

0 

0 

H^Ui(an) 

0 

0 


(19) 

Now let us find the inverse matrix (Pfc) ^ of the matrix P^. If we denote ith 
row of the inverse matrix {Pk)~^ by {Pk)^^, then we have 


{Pk)i = 


(k-i) (k-i)' 

YlrUo{otr) 0...0 ri^ii~^Ui{ar) ... n-k jar) 0... Oj 

( 20 ) 

i = l,l + fc,l + 2fc,...,n — fc + l;r = ^ 


{Pk)^ = 


(fc-l) (fc-2) 

0 rj^Uoiar) 0...0 ... U n-k (ar) 0...0 

1 . „ _ i+k—2 


-iT 


( 21 ) 


i = 2, 2 + fc, 2 + 2fc,..., n — fc + 2; r = ^ 


(fc-l) (fc-l) (fc-3) 

(l^fe)i=lo 0 T]rUo{ar) 0^^ ... o”^ o”^ 

z = 3, 3 + fc, 3 + 2fc,..., n - fc + 3; r = 


( 22 ) 


iPk)i 



(fc-l) 

r]^Uo{ar) 0...0 


(fc-l) 

0...0 n„u~^^U n-k (ar) 


T 


Oj ’ 

(23) 
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i = A: — 1,2A:—1,3A:—l;r = 


' (^-1) (fc-i) (^-1) 1 

0^ ri^Uo{ar) 0^ ... 0^ (a^) ’ 

i = k, 2k, 3k,... ,n;r = ^ 

where /k = |, and for r = 1, Thus, we obtain 


(Pk)-^ 


r]iUo{ai) 

0 0 

rj^fx iUi{ai) ••• 

0 

rj^Uoiai) 0 

0 

0 

0 0 

?7iC/o(ai) 

V2Uo{a2) 

0 0 

r72Ai“5[/i(a2) ••• 

0 

T]2Uo{a2) 0 

0 

0 

0 0 

%C^o(a2) 

r7nC/o(an) 

0 0 

■■• ?7n/x-5[/i(an) ... 

0 

r]nUoia^) 0 

0 

0 

0 0 

VfUoia^) 






0 


0 

n — k 


0 0 
0 q2t^~^^~^U,i-k _^{a2) 


2jrUn^{ai) 

0 

0 


q2^l ^ ^Ur^_^{a2) 


qnq, 


0 

n — k 

?7„/^-Tir[7^(an) 


0 0 ••• (tt") 

It ' ' fc 

(25) 

By combining (ITHl) . (fT^ and (l25|) and using the equality {AkY = Pk{JkY{Pk) ^ 
|12] . we compute the Zth powers of the matrix as 


[AkY = Pk[JkY{Pk)-^ =Q[l) = [q^J (0) • 


(26) 


So, for i,J = l,n 



1 0, 

if other, 


9y(0-j 

[ar)Uj-vij (or)) if i = j and i 

= j + k,j- 

V r=l fc 

k 

(27) 

r k 

; i,j = 0 mod [k) , 


II 

; i,j = m mod (fc) , 

(28) 


where m = 1, fc — 1. 

Corollary 2 Let A^ he n—square [n = 2ks, s G N; a, 6, c G R \ {0}) 
[2k + 1) —diagonal Toeplitz matrix as in flj), from Theorem 1 

a^2Vbccos(^-^^ (29) 

(r = 1, ^). In that case, there exists the inverse and negative integer powers of 
the matrix A^. 
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4 Numerical Examples 

Example 3 Setting n = 8, for k = 1 and I = 8, we have 

Ji = (izag(Ai, A2, A3, A4, As, Ae, A7, As) 

= diag{a — 1.879-\/^, a — 1.532'/bc, a — Vbc, a — O.MlVbc, 
a + 0.347-\/fc, a + Vbc, a + 1.532-\/^, a + 1.879'/bc) 

and 

[Aif = Pi(Ji)3(Pi)-i = Q(3) = 

+ 8abc 8afb + 2&^c 3a&^ b^ 

8afc + 2bc^ + 6abc 3a^b + 3b^c 3ab^ 
3ac^ 3a^c + 3bc^ + 6abc 3a^b + 3b^c 
3ac^ 3afc + 3hc^ + 6abc 
(gij(3)) =0 3ac^ 3a^c + 3bc^ 

0 0 c3 3ac^ 

0 0 0 c3 

0 0 0 0 


0 0 0 0 

53 0 0 0 

3ab^ 53 0 0 

3a^b + 3b^c 3ab^ b^ 0 

+ 6abc 3afb + 36^c 3ab’^ b^ 

3afc + 3bc^ + 6a6c 3afb + 35^c 3a6^ 

3ac^ 3afc + 3bc^ + 6abc 3afb + 25^c 
c3 3ac^ 3afc + 2bc^ a^ + 3abc 

For k = 2J = 4:,a = i — l,b = i + A and c = i — 2 

J2 = (izo(/(Ai, Ai, A2, A2, A3, A3, A4, A4) 

= diag{-lA3& - 3.884i, -1.536 - 3.884i, -1.205 - 0.865i, -1.205 - 0.865i, 
-0.795 + 2.865i, -0.795 + 2.865i, -0.464 + 5.884i, -0.464 + 5.884i) 

and 

(A2)^ = P2(J2)"(P2)-' = Q(0 = 

174 + 361 0 336- 256Z 0 

0 174 + 36i 0 336- 256i 

-48 + 224i 0 429 + 36i 0 

0 -48 + 224i 0 429 + 36z 

-105 + 90Z 0 -60 + 340Z 0 

0 -105 + 90Z 0 -60 + 340i 

-36-52i 0 -105 + 90Z 0 

0 -36- 52i 0 -105 + 90Z 
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-357 - 306i 0 -396 + 20* 

0 -357- 306* 0 

492 - 404i 0 -357 - 306* 

0 492 - 404i 0 

429 + 36* 0 336 - 256* 

0 429 + 36* 0 

-48 + 224i 0 174 + 36* 

0 -48 + 224i 0 


0 

-396 + 20i 
0 

-357 - 306* 
0 

336 - 256* 
0 

174 + 36* 


For k = A 

■h 


diag{Xi, Ai, Ai, Ai, A 2 A 2 , A 2 , A 2 ) 


where Xi = a — Vbc, X2 = a + Vbc and 

(AiY = PiidiYiPi)-^ = Q{i) = 


qiiil) 912(0 913(0 914(0 915(0 9i6(0 917(0 9i8(0 

92l(0 922(0 923(0 924(0 925(0 926(0 927(0 928(0 

93l(0 932(0 933(0 934(0 935(0 936(0 937(0 938(0 

/ /.NN _ 941 (0 942(0 943(0 944(0 945(0 946(0 947(0 948(0 

95l(0 952(0 953(0 954(0 955(0 956(0 957(0 958(0 

96l(0 962(0 963(0 964(0 965(0 966(0 967(0 968(0 

97l(0 972(0 973(0 974(0 975(0 976(0 917(0 978 (0 

98i(0 982(0 983(0 984(0 985(0 986(0 987(0 988(0 

9ii(0 = 922(0 = 933(0 = 944(0 = 955(0 = 966(0 
= 977(0 = 988(0 = 0.5(a + y/bcY + 0.5(a — 

915(0 = 926(0 = 937(0 = 948(0 = [0.5(a + VtcY - 0.5(0 - 

95i( 0 = 962(0 = 973(0 = 984(0 = [0.5(a + - 0.5(o- VbcYjfx^ 

912(0 = 913(0 = 914(0 = 916(0 = 917(0 = 918(0 = 0 

921(0 = 923(0 = 924(0 = 925(0 = 927(0 = 928(0 = 0 

93l(0 = 932(0 = 934(0 = 935(0 = 936(0 = 938(0 = 0 

941 (0 = 942(0 = 943(0 = 945(0 = 946(0 = 947(0 = 0 

952(0 = 953(0 = 954(0 = 956(0 = 957(0 = 958(0 = 0 

96i(0 = 963(0 = 964(0 = 965(0 = 967(0 = 968(0 = 0 

97i(0 = 972(0 = 974(0 = 975(0 = 976(0 = 978(0 = 0 

98i(0 = 982(0 = 983(0 = 985(0 = 986(0 = 987(0 = 0 


where /* = f • 

Example 4 Setting n = 6 , for k = 1, I = —3, a = i, b = i + 1 and c = i — 
Ji = diag{Xi, X 2 , Xs, Xi, X 5 , Xq) 

= diag{-1.5A8i, -0.763*, 0.371*, 1.629*, 2.763*, 3.548*) 
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and 




i^iji 3)) 


= Pi(Ji)-3(Pi)-i =Q(-3) = 

5i -1.286- 1.286i 

1.286 - 1.286Z 0.429i 

4.571 -1.571+ 1.57H 

-2.857- 2.857i 1.714 

-2.857i 1.143+ 1.143i 

-4 + 4i -2.857i 


-4.571 

1.571 + 1.57H 
3.286i 

2.429 - 2.429Z 
1.714 

-2.857- 2.857i 


2.857- 2.857i 
1.714 

-2.429- 2.429i 
3.286i 

-1.571 + 1.571i 

4.571 


-2.857i 

-1.143+ 1.143Z 
-1.714 

1.571 + 1.571i 
0.429* 

1.286- 1.286 


4 + 4* 
-2.857i 
2.857- 2.857* 
-4.571 

-1.286- 1.286* 
5* 


For k = 2, I = —4, a = 3, 6=1 and c = 6 


and 


Ji = d*ag(Ai, Ai, A 2 , A 2 , As, As) 

= diag(-0.464, -0.464, 3,3,6.464,6.464) 



= P2{J2)-\P2) 

1 

II 

1 

4) = 




5.395 

0 

-3.112 




0 

5.395 

0 



= 

-18.667 

0 

0 

-18.667 

10.778 

0 




32.296 

0 

-18.667 




0 

32.296 

0 





0 

0.897 

0 




-3.112 

0 

0.897 




0 

-3.112 

0 




10.778 

0 

-3.112 




0 

5.395 

0 




-18.667 

0 

5.395 


5 Complex Factorization 

The well-known Fibonacci polynomials F(x) = are defined in 

by the recurrence relation, 


Fn(x) = xFn-i(x) + Fn- 2 {x) (** > 3) 
where Fo{x) = 0 , Fi{x) = I and ^ 2 ( 0 ;) = x. 


12 






Corollary 5 Let the matrix Ak be n—square (n = 2ks, s G N) k—diagonal 
Toeplitz matrix as in m with a := x,b := i and c := i where i = \/—1. Then 

det(Afe) = (Fn+i(a;))'= (30) 

where (x) are n—th Fibonacci polynomial. 
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